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Abstract 

A geometric setup for control theory is presented. The argument is devel- 
oped through the study of the extremals of action functionals defined on 
piecewise differentiable curves, in the presence of differentiablc non-holonomic 
constraints. Special emphasis is put on the tensorial aspects of the theory. 
To start with, the kinematical foundations, culminating in the so called vari- 
ational equation, are put on geometrical grounds, via the introduction of the 
concept of infinitesimal control. On the same basis, the usual classification of 
the extremals of a variational problem into normal and abnormal ones is also 
rationalized, showing the existence of a purely kinematical algorithm assigning 
to each admissible curve a corresponding abnormality index, defined in terms 
of a suitable linear map. The whole machinery is then applied to constrained 
variational calculus. The argument provides an interesting revisitation of 
Pontryagin maximum principle and of the Erdmann-Weierstrass corner con- 
ditions, as well as a proof of the classical Lagrange multipliers method and 
a local interpretation of Pontryagin's equations as dynamical equations for a 
free (singular) Hamiltonian system. As a final, highly non-trivial topic, a suf- 
ficient condition for the existence of finite deformations with fixed endpoints 
is explicitly stated and proved. 
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1 Introduction 



Constrained variational calculus has been extensively studied since the beginning 
of the XX century and has been recently revived by its applications to control 
theory. Among others, we mention here the pioneering works of Bolza and Bliss [T], 
the contribution of Pontryagin [Ij and the more recent developments by Sussman, 
Agrachev, Hsu, Montgomery and Griffiths [18 1 121 ^ 12 ^ 120 1 112j. characterized by the 
use of the mathematical instruments provided by differential geometry. 

As it is already clear in Bliss work, the primary task is finding out to what 
extent the classical tools for the determination of constrained extremals in ordinary 
Calculus (e.g. Lagrange multipliers) may be extended to the study of functionals. 

Quite naturally, the argument requires a preliminary characterization of the 
admissible curves and of their deformations. In this connection, the analysis reveals 
the existence of curves having a pathological behaviour with respect the fixed-end 
deformations, an aspect well known in the literature, under the name of abnormal 
extremals [Ml^ . 

In this paper we propose a fresh approach to the argument. We shall deal with 
systems described by a finite number of variables g^, . . . subject to a set of 
differentiable conditions of the form 

^ = ^P\t,q\...,q\z\...,zn (1.1) 

expressing the derivatives ^ in terms of a smaller number of control variables 
z^, A = 1, . . . ,r. Every set of functions = g*(t), = z'^{t) consistent with 
the requirement (jl.ip will be called an admissible evolution of the system. 

In addition to the constraints p.2|) we shall also consider an action functional 

1= r L{t,q\...,q^,z\...,z'^)dt (1.2) 
Jto 

expressed as the integral of a suitable "cost function", or Lagrangian L(t,q,z) 
along the admissible evolutions of the system. Within the stated context, we shall 
discuss the (local) extremals of I with respect to the admissible deformations 
leaving the endpoints q^{to), g*(ti) fixed. 

Geometrically, an intrinsic formulation of the problem is obtained introducing 
a fibre bundle Vn+i with coordinates t,q^ , . . . ,q^ ^ called the event space, 

denoting by ji(V„+i) ^ V^+i the associated first jet bundle, and regarding the 

control equations (jl.ip as the representation of a submanifold A ^ ii(Vn+i) 
describing the totality of admissible kinetic states of the system. In the resulting 
environment, a control for the system is then simply a section a : Vn+i A , 
locally represented as z"^ = z"^(t , g-*^, . . . , g") . 

The infinitesimal deformations of an admissible section 7 : M ^ Vn+i are 
discussed via a revisitation of the familiar variational equation. The novelty of the 
approach relies on the introduction of a transport law for vertical vector fields along 
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7, yielding a covariant characterization of the "true" degrees of freedom involved 
in the description of the most general admissible infinitesimal deformation. 

The analysis is subsequently extended to arbitrary piecewise differentiable evo- 
lutions 7 consisting of families of contiguous closed arcs 7^*^ : [0^-1,0^] — >■ Vn+i, 
to = ao < ai < ■ ■ ■ < aj\f = ti denoting any (finite) partition of the interval [to,ti]- 

In this connection, special attention is paid to the infinitesimal deformations 
vanishing at the endpoints. No restrictions are posed on the deformability of the 
intervals [a^-i, as] or on the mobility of the "corners" 7(0^) ,s = l,...,A'" — 1. 

The argument allows to assign to every admissible evolution a corresponding 
abnormality index, rephrasing and putting on geometrical grounds the traditional 
attributes of normality and abnormality commonly found in the literature [181 [20] . 

Equally important, in a variational context, is finding out which infinitesimal 
deformations are actually tangent to finite deformations with fixed end points. 
To account for this aspect, the admissible evolutions will be classified into ordinary 
ones, if every admissible infinitesimal deformation vanishing at the endpoints is 
tangent to some finite deformations with fixed end points, and exceptional ones 
in the opposite case. An important result, proved in Appendix B and extending 
a result proved in [2^ in the case of linear constraints, is the relationship be- 
tween abnormality index and ordinariness, resulting in the fact that every normal 
evolutions is also, automatically, an ordinary one. 

After these preliminaries, attention is focussed on the study of the extremals of 
the functional (II. 2p . Once again, the approach relies on a fully covariant algorithm, 
summarizing the content of Pontryagin's maximum principle 01 [19] and of the 
Erdmann Weierstrass corner conditions [H [T^ . The resulting equations are shown 
to provide sufficient conditions for any evolution, and necessary and sufficient 
conditions for an ordinary evolution 7 to be an extremal. The same setup is seen 
to provide an elegant and concise proof of the Lagrange multipliers method. 

In the final part of the paper, the geometric content of the algorithm is further 
enhanced, lifting everything to a fiber bundle C{A) — > A, here called the contact 
bundle, defined as the pull-back of the phase space V*{Vn+i) through the fibered 
morphism 

C{A) . V*{Vn+l) 

A > Vn+l 

In the resulting framework, the ordinary extremals of the original variational 
problem are seen to arise as projections of the solutions of a corresponding free 
variational problem on C{A). The Euler-Lagrange equations for the newer problem 
are written in terms of the so called Pontryagin Hamiltonian TC := pi ip^'—L, viewed 
as a function on C{A). The same environment provides also a characterization of 
the ordinary abnormal evolutions of the system as projections of the extremals of 
a purely geometric action principle based on the Liouville 1-form of C{A). 

The argument is completed by a discussion of the meaning of the Pontryagin 
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Hamiltonian. It is seen that, under suitable regularity assumptions, the original 
constrained variational problem is locally equivalent to a free Hamiltonian problem 
on the phase space. 



2 Geometric setup 
2.1 Preliminaries 

In this Section we outline the construction of a geometrical setup especially suited 
to the development of constrained variational calculus and smooth control theory. 
In this connection, see also [Il[2l[ll[5l[7l[8l[9l[l2l[l8l[l9l[2T]and references therein. 

(i) Let V„+i ^ ffi denote a fibre bundle over the real line, henceforth called the 
event space, and referred to local fibered coordinates t, g^, . . . , g". 

Every section 7 : M — > Vn+i , locally described as = q^{t) , will be interpreted 
as an evolution of an abstract system *B, parameterized in terms of the independent 
variable t. The terminology is borrowed from mechanics, where !B is identified with 
a material system, Vn+i with the configuration space-time of 55 , and t : Vn+i — > K 
with the absolute time function. 

Pursuing this analogy, the first jet bundle ji{Vn+i), referred to local jet coor- 
dinates t,q'^,q'^, will be called the velocity space. 

As reported in Appendix A, the projection ji(Vn+i) — ^ Vn+i makes ji{Vn+i) 
into an affine bundle, modelled on the vertical space V{Vn+i)- The vertical bundle 
associated with this projection will be denoted by l^(ji(V„,+i)) . 

Every section 7 : M ^ V„+i admits a lift ji(7) : M ji(Vn+i), expressed 
in coordinates as = Q^i^) j = ^ ■ In a similar way, every vertical vector 
field X = over V„+i may be lifted to a field J(X) = X'^ + X'^ over 

ji(V„+i), with X* := + q'^^F- Both arguments are entirely standard (see 
e.g. |16j )■ and will be regarded as known. 

(ii) The restriction of the vertical space V{Vn+i) to the section 7 determines a 
vector bundle V{'y) M, called the vertical bundle over 7. 

Proposition 2.1 The first jet space ji{V{^)) is canonically isomorphic to the 
vector bundle over R formed by the totality of vectors Z along ^1(7) satisfy- 
ing (^Z,dt) = 0. With this identification, the fibration vr* : ji{V{'y)) ^(7) 
coincides with the restriction to Ji(^(7)) of the push-forward of the projection 

TT : ji(Vn+l) Vn+l- 

Proof. Given t* E M and a section X : M — > V{'j) , choose any vector field X 
defined in a neig hborhood U B -f{t*) and satisfying = X{t) Vt e 7"^^)- 

In coordinates, setting ^ : q" = g*(t) , X = X^{t) (^)^, the lift of the field X 
at the point ji{'j){t*) takes the value 



t=t' /ii(7){i-) 
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All assertions of Proposition 12.11 follow immediately from this fact. □ 

Consistently with eq. (12. ip . given any section X : R ^ ^(7) ! the jet extension 
ji {X) will be called the lift of X to the curve ji (7) . In local coordinates, eq. (|2.1|) 
provides the representation 



n{x) = XM — 1 + 



d_ 



dX' / d 



(2.2) 



Every local coordinate system i, in Vn+i determines fibered coordinates t, 
in V{^) and t,u^,u^ in ji{V{'y)), based on the identifications 



X = u\X) 



_d_ 

d 



hiimz)) 



V X G 1/(7) (2.3a) 
VZGii(y(7)) (2.3b) 



In terms of these, the jet-extension of a section = n*(t) takes the standard form 
= u^{t) , ii^ = ^ , while the projection vr* : jiiVin)) ^(7) is described by 
■u*(7r*(Z)) = n*(Z) . In particular, the vertical subbundle V{ji{'y)) coincides with 
the submanifold of ji{V{'y)) locally described by the equation li* = 0, i = 1, . . . ,n. 



Corollary 2.1 The vector bundles V{ji{'j)) 
isomorphic 



and 1^(7) — ^ M are canonically 



Proof. As pointed out in Appendix A, for each z G ji(Vn+i) the affine character 
of the fibration ji(Vn+i) — ^ Vn+i determines an isomorphism between the vertical 
spaces Vz{ji{Vn+i)) and T4-(2)("^'n+i) j expressed in coordinates as 



_d_ 



d_ 



7r(z) 



In particular, for z = ii(7)(t) , our previous definitions imply the identifications 
7r(z) = 7(t) , K(.)(V„+i) = F(7)|i , V,{h{Vn+i)) = V{h{l))\t ■ By varying t, this 
gives rise to a vector bundle isomorphism 



V{ji{l)) 
t 



y(7) 
t 



(2.4) 



expressed in coordinates as 



(2.5) 
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□ 



(iii) Together with the vertical bundle ^^(7) — > M it is worth considering the dual 
bundle V*{j) R, identical to the pull-back on 7 of the phase space V*{Vn+i) ■ 
With the notation of Appendix A, the situation is expressed by the commutative 
diagram 

V*{l) > V*{Vn+l) 

t TT (2.6) 

M Vn+1 

The elements of V*{'~f) will be called the virtual l-forms along 7. 

Notice that, according to the stated definition, a virtual 1-form A at a point 
7(t) is not a 1-form in the ordinary sense, but an equivalence class of 1-forms 
under the relation 

A ~ A' ^ A - A' oc (dt)^(t) (2.7) 

For simplicity, we preserve the notation ( , ) for the pairing between V{^) and 
V*{'y). Also, given any local coordinate system t,q^ in V„+i , we refer V*{'y) to 
fiber coordinates t,pi , with Pi{X) = (A, • The virtual 1-forms along 

7 determined by the differentials dg* will be denoted by i = 1, . . . , n. 

For completeness we remark that, according to eq. (I2.6p . each fiber V*{^)\t 
is isomorphic to the subspace of the cotangent space T*^^s^{Vn+i) annihilating the 
tangent vector to the curve 7 at the point 7(t). Formally, this viewpoint is im- 
plemented by setting u)* = [dq^ ~ ^ ■ Although apparently simpler, this 
characterization of V*{'j) has some drawbacks in the case of piecewise differen- 
tiable sections. We shall therefore stick to the original definition. 

(iv) Let us now see how the geometric setup developed so far gets modified in 
the presence of non-holonomic constraints reducing the number of independent 
velocities. Under the explicit assumption of differentiahle constraints, the situation 
is summarized into a commutative diagram of the form 

A Jl(Vn+l) 

TT TT (2-8) 

Vn+1 ^=^= Vn+1 

where 

• A ^ Vn+1 is a fiber bundle, representing the totality of admissible velocities] 

• the map A ^ ji(Vn+i) is an imbedding; 

• a section 7 : M — > Vn+i is admissible if and only if its lift ^1(7) factors 
through A, i.e. if and only if there exists a section 7 : M — > ^ satisfying 

hil) = i • 7- 
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The geometry of the submanifold A is weh known from non-holonomic me- 
chanics. Concepts hke the vertical bundle V{A) and the contact bundle C{A) will 
be freely used in the following. Their properties are briefly reviewed in Appendix A. 
For a detailed analysis, see e.g. [22j and references therein. 

Every section 7 : M — s- ^ satisfying i • 7 = ji(vr • 7) will be called admissible. 
With this terminology, the admissible sections of A are in 1-1 correspondence 
with the admissible evolutions 7 : M ^ Vn+i . The manifold A will be referred to 
local fibered coordinates t, q^, . . . , g", z\ . . . , z*", with transformation laws 

t = t, r=nt,q\...,qn, z^ = z^{t,q\...,q^,z\...,zn (2.9) 

The imbedding i : A^ jiO^n+i) will be locally expressed as 

q'=^P\t,q\...,q-,z\...,zn i = l,...,n (2.10) 



with rank 



■■■■4,") 

<9(2l--2'-) 



= r. With this notation, given any section 7 : M ^ ^ 

described in coordinates as = (?*(t) , z^ = z^{t) , the admissibility requirement 
i • 7 = ji (vr • 7) takes the explicit form 

^ = r{t,q\t),...,q^{t),z\t),...,z^{t)) (2.11) 

Eqs. (|2.11|) indicates that, for any admissible evolution of the system, the 
knowledge of the functions z^{t) determines q^{t) up to initial data. 

On the other hand, in the absence of specific assumptions on the nature of 
the manifold A, such as e.g. the existence of a canonical factorization of the form 
-2^ Xr Vn+i , 2 denoting an (r + l)-dimensional fiber bundle over M with coordi- 
nates t,z^, . . . ,z^ , the functions z'^{t) , in themselves, have no invariant geometrical 
meaning. 

To pursue the idea of the z^^s as "handles" controlling the evolution of the 
system, attention should rather be shifted on sections a : Vn+i A . 

Every such section will be called a control for the system; the composite map 
i ■ a : Vn+i ii(Vn+i) will be called an admissible velocity field. 

In local coordinates we have the representations 



z 



A 



I ■ a 



z^(t,q\...,q^) (2.12a) 
r{t,q\...,q'',z^{t,q\...,q^)) (2.12b) 



confirming that the knowledge of a does indeed determine the evolution of the 
system from any given initial event x € V„+i through a well posed Cauchy problem. 

A section 7 : M ^ Vn+i and a control a : Vn+i A will be said to belong to 
each other if and only if the lift 7 : M — > ^ factors into 7 = a • 7, i.e. if and only if 
the jet extension ji{'y) coincides with the composite map i • u • 7 : M — > ji(Vn+i). 

(v) Given any admissible section 7, let ^(7) — > M denote the vector bundle 
formed by the totality of vectors along 7 annihilating the 1-form dt. For each 
t G M, the fiber of A (7) over t will be denoted by ^(7)1^ . 
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On account of Proposition 12. H the push-forward : T{A) T{ji(Vn+i)) 
gives rise to a bundle morphism 

(2.13) 

vh) — vii) 

making ^(7) into a subbundle of Ji(V^(7)), fibered over ^^(7). 

Once again all arrows in diagram (|2.13p . regarded as maps between vector 
bundles over M, have the nature of homomorphisms. The kernel of the projection 
^(7) ^(7); clearly identical to the restriction of the vertical bundle V{A) 
to the curve 7, will be denoted by ^(7), and will be called the vertical bundle 
along 7 . 

Every fibered coordinate system t,q^,z^ in A induces coordinates t,u^,v^ in 
A (7) according to the prescription 

^ = "'^^^(^) -+"^(^^(9^) - Vie ^(7) (2.14) 

In terms of these, and of the jet coordinates t, u*, on Ji(l^(7)) , the morphism 
(|2.13p is locally described by the system 

d'4^'\ k , ( dV\ A 



t = t, u^ = u\ '''=[j^^)f (2.15) 

while the vertical bundle ^(7) coincides with the slice = in A (7). 

(vi) For later use, let us finally examine the "constrained" counterpart of diagram 
(|2.4p . To this end we observe that the morphism (I2.13P maps V (7) into the vertical 
subbundle ^(ji(7)) C Ji(^(7)) • Composing with the morphism (j2.4p . and setting 
Q ■= Q ■ i* , this gives rise to an injective homomorphism 

(2.16) 



In coordinates, eqs. ()2.5p . ()2.15p provide the representation 



d 



_d_ 



(2.17) 
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2.2 Deformations 



Quite generally, given a fiber bundle M — > N and a section u : N —)■ M, a 
deformation of is a 1-parameter family of sections z^^ , ^ G {—e,e) depending 
differentiably on ^ and satisfying vq = v. 

For each y (z N , the curve ^ — > J^g(y) is called the orbit of the deformation 
through the point v^y) . The vector field along v tangent to the orbits at ^ = is 
called the infinitesimal deformation associated with . 

Let us apply all this to the case in study. 

(i) According to the stated definition, a deformation of a control a : Vn+i ^ ^ is 
a one parameter family of controls depending differentiably on ^ and satisfying 
(To = (T. In coordinates, preserving the representation (12. 121 a). we shall write 



z 



^-C^{i,t,q\...,q^) (2.18) 



with ('"^(O, i , g^, . . . , g") = z'^^t, q^, . . . , g") . The infinitesimal deformation associ- 
ated with cTg is expressed locally as (af^)cr ' -^'^o™ this it is readily seen 
that the totality of infinitesimal deformations of a control a coincides with the 
totality of vertical vector fields along a . 

(ii) When deforming an admissible section 7 : M ^ Vn+i care must be taken of the 
presence of constraints. Given any such 7, a deformation 7^ is called admissible 
if and only if each section 7^ : M ^ V^+i is admissible. In a similar way, a 
deformation 7^ of an admissible section 7 : M ^ ^ is called admissible if and only 
if all sections 7^ : M — > ^ are admissible. 

As pointed out in § 2.1, the admissible sections 7 : M — > V„+i are in 1-1 
correspondence with the admissible sections 7 : M — > ^ through the relations 

7 = 7r-7, ji(7) = i.7 (2.19) 

Every admissible deformation of 7 may therefore be expressed as 

7^ = vr • 7g 

7^ : M ^ ^ denoting an admissible deformation of 7. 

In coordinates, preserving the representation 7:9* = q^{t) , z^ = z^{t) , the 
admissible deformations of 7 are described by equations of the form 

ir- q' = v'{i,t), z^ = C^{^,t) (2.20) 

subject to the conditions 

^\0,t)=q^{t), C^{0,t) = z^{t) (2.21a) 

^ = ij\t,^\(,t),CH(,t)) (2.21b) 
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Setting 



the infinitesimal deformation tangent to 7^ is described by the vector field 

^ = X'it)(4j] (2.22) 

7 



while eq. ()2.21l b) is reflected into the relation 



dt dtdi 



5=0 



^ + hrr ^ (2-23) 



commonly referred to as the variational equation. The infinitesimal deformation 
tangent to the projection 7^ = tt • 7^ is similarly described by the field 

Collecting all previous results and recalling the definitions of the vector bundles 
1^(7) and A (7) we get the following 

Proposition 2.2 Let 7 : M ^ Vn+i o-nd 7 : M ^ ^ denote two admissible 
sections, related by eq. 112. 19\) . Then: 

i) the infinitesimal deformations of 7 and of 7 are respectively expressed as 
sections X : M ^ 1/(7) anrf X : M ^ ^(7) ; 

ii) a section X : M ^ ^(7) represents an admissible infinitesimal deforma- 
tion of 7 if and only if its first jet extension factors through ^(7), i.e. if 
and only if there exists a section X : M — > ^(7) satisfying ji{X) = i^^X; 
conversely, a section X : M — > ^(7) represents an admissible infinitesimal 
deformation of 7 if and only if it projects into an admissible infinitesimal 
deformation of 7, i.e. if and only if i^,X = ji['K^X) . 

The proof is entirely straightforward, and is left to the reader. 

From a structural viewpoint, Proposition 12.21 establishes a complete symme- 
try between the roles of diagram ()2.8p in the study of the admissible evolutions 
and of diagram (j2.13p in the study of the admissible infinitesimal deformations, 
thus enforcing the intuitive idea that the latter context is essentially a "linearized 
counterpart" of the former one. 



2.3 Infinitesimal controls 

According to Proposition 12.21 the admissible infinitesimal deformations of an ad- 
missible section 7 : M — > Vn+i are in 1-1 correspondence with the sections 
X : M ^ ^(7) satisfying the consistency requirement i^^X = ji(7r*X) . 
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In local coordinates, setting X = X^{t) ^ +T'^{t) > the stated requirement 
is expressed by the variational equation 

all coefficients being evaluated along the curve 7 . 

Exactly as it happened with eq. (12. lip , eq. ()2.25p indicates that, for each 
admissible X, the knowledge of the functions T^{t) determines the remaining 
X'^(t) up to initial data, through the solution of a well posed Cauchy problem. 

Once again, however, the drawback is that the components F^, in themselves, 
have no invariant geometrical meaning, but obey the non-homogeneous transfor- 
mation law 

under arbitrary coordinate transformation. The difficulty is overcome introducing 
a linearized version of the idea of control. 

Definition 2.1 Let 7 : M — > Vn+i denote an admissible evolution. Then: 

• a linear section h : ^^(7) ^(7); meant as a vector bundle homomorphism 
satisfying tt^ ■ h = id, is called an infinitesimal control along 7; 

• the image 'H{^) := h{V{'y)), viewed as a vector subbundle of ^(7) M, 
is called the horizontal distribution along 7 induced by h; every section 
X : M — > ^4(7) satisfying X(t) € 7^(7) V t G M is called a horizontal section. 

Remark 2.1 The term infinitesimal control is intuitively clear: given an admis- 
sible section 7, let a : Vn+i — ^ A denote any control belonging to 7, i.e. satisfying 
£7 • 7 = 7. Then, on account of the identity vr,,, • cr* = (vr • a)* = id, the tan- 
gent map : T{Vn+i) T{A), restricted to V{'y), determines a linear section 
(T* : ^(7) ^(7)- The infinitesimal controls may therefore be thought of as 
equivalence classes of ordinary controls having a first order contact along 7. 



Given an infinitesimal control h : V^j) — > ^(7), on account of Definition 12.11 
it is easily seen that the horizontal distribution "^(7) and the vertical subbundle 
^(7) split the vector bundle ^4(7) into the fibered direct sum 

A (7) = 7^(7)01,^(7) (2.26) 

thus giving rise to a couple of homomorphisms Vh '■ ^(7) — *■ ^^(7) (horizontal 
projection) and Vv : ^(7) — *■ ^(7) (vertical projection), uniquely defined by the 
relations 

Vh = h-TT^ ; Vv = id-Vu (2.27) 

In fiber coordinates, preserving the notation ()2.3l a). (I2.14p . every infinitesimal 
control h : V{'j) ^(7) is represented by a linear system of the form 

v^ = hi'^{t)u' (2.28) 
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In this way: 

• the horizontal distribution 7^(7) is locally spanned by the vector fields 



di := h 



_d_ 



_d_ 



+ hi 



_d_ 



(2.29) 



every vertical vector field X = X*(t)(^) along 7 may be lifted to a 
horizontal field h{X) along 7, expressed in components as 



h{X) = X\t)di = X\t) 



d_ 

dq' 



d 



(2.30) 



every vector X = X^ 



' a \ 



G ^(7) admits a unique represen- 



tation of the form X = Vh{X) + Vv{X) , with 



VH{X)=X'd,, VviX) 



X'h,, 



d 



_d_ 

dz^ 



(2.31) 



The role of Definition 12.11 in the study of the variational equation (I2.25|) is 
further enhanced by the following 

Definition 2.2 Let h be an infinitesimal control for the (admissible) section 7. 
A section X : M — > ^^(7) is said to be h-transported along 7 if and only if its 
horizontal lift h{X) : M ^(7) is an admissible infinitesimal deformation of j , 
i.e. if and only if • h{X) = ji{X) . 

In coordinates, setting X = X^{t) (^)^ and recalhng eqs. (I2.25|) . (I2.30p . the 
condition for /i-transport is expressed by the linear system of ordinary differential 
equations 



dX' 
IT 



dq'' 



+ hk 



X' 



(2.32) 



From the latter, recalling Cauchy theorem, we conclude that the /i-transported 
sections of 1^(7) form an n-dimensional vector space Vh, isomorphic to each fibre 
y(7)|( through the evaluation map X ^ X{t) . We have thus proved: 

Proposition 2.3 Every infinitesimal control h : V{'^) ^(7) determines a triv- 
ialization of the vector bundle V^j) ^ M. 



Proposition 12.31 provides an identification between sections X : R ^ ^(7) and 
vector valued functions X : R — > V/i , and therefore, by duality, also an identifica- 
tion between sections A : R — > y*{l) and vector valued functions 
A : R — > y^*, thus allowing the introduction of an absolute time derivative 
for vertical vector fields and virtual 1-forms along 7. 

The algorithm is readily implemented in components. To this end, let {e(a)}, 

{e^"')} denote any pair of dual bases for the spaces Vh, V^. By definition, each e(a) 
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is a vertical vector field along 7, obeying the transport law (j2.32p . 

'(a) \dqij-y 

de 



In coordinates, setting e^^-j = e,\ (^) , this implies the relation 



^ = %^)9.V'^ (2.33a) 

In a similar way, each e^"^ is a contact 1-form along 7, expressed on the basis cj* 
as e(") = e^f Qj\ with ef^ e(^) = . 

On account of eq. (j2.33l a). the components e'f' obey the transport law 



The functions 

. , , de,\ 

- dt - dt ^^-^^^^ 

will be called the temporal connection coefficients associated with the infinitesimal 
control h in the coordinate system t,^*. Comparison with eqs. ()2.29p . ()2.33l a.b) 
provides the representation 

Given any vertical vector field X = along 7, the definition of the 

operator is summarized into the expression 

d / V {a)\ d I (a) .A d ( . (a)\ j ( d 

written more simply as 

DX ( dX^ ^ .\( Q \ , . 

+ XV,n — (2.35a) 



Bt \ dt ) \dqj 

with the coefficients Tj-' given by eq. (|2.34l a). In a similar way, given any virtual 
1-form A = Aj u)* , the same argument provides the evaluation 



After these preliminaries, let us go back to the variational equation (j2.25p . By 



means of the projections p.3ip . every section X = X*(^). +T^[^^). splits 



into the sum 

X = Vh{X) + Vv{X) = h{X) + Y (2.36) 
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w\ihX = TT,{X),Y = Vv{X) = {T^-hi^X'){^).. 

On the other hand, on account of eq. (I2.29p . the variational equation ()2.25p is 
mathematicaUy equivalent to the relation 



dt 

Recalling eqs. (I2.34l b). ()2.35l a). ()2.36p . as well as the representation (12.170 of the 
homomorphism V{^) ^(7)) the latter may be written synthetically as 

DX 



q[Y) = q{Vv{X)) (2.37a) 
A) 



Dt 

or also, setting X = X"e(a)) Y = Y"^ (^-^x) and expressing everything in com- 
ponents in the basis e(a) 

dX"" 



dt 



(e('^),^(y))=e(r)(|5) (2.37b) 



Exactly as the original equation (|2.25p . eq. (|2.37l a) points out that every in- 
finitesimal deformation X is determined by the knowledge of a vertical vector field 
Y = Y^ (af^').^. through the solution of a well posed Cauchy problem. 

The advantage is that, in the newer formulation, all quantities have a precise 
geometrical meaning relative to the horizontal distribution 7^(7) induced by the 
infinitesimal control h. 

On the other hand, one cannot overlook the fact that, in the standard formula- 
tion of control theory, no distinguished section h : 1^(7) — > ^(7) is provided, and 
none is needed in order to formulate the results. In this respect, the infinitesimal 
control h plays the role of a gauge field, useful for covariance purposes, but un- 
affecting the evaluation of the extremals. Accordingly, in the subsequent analysis 
we shall employ /i as a user-defined object, eventually checking the invariance of 
the results under arbitrary changes /i — > /i'. 



2.4 Corners 

To complete our geometrical setup we have still to consider the fact that, in general, 
control theory does not deal with sections in the ordinary sense, but with piecewise 
differentiahle evolutions, defined on closed intervals. To account for this aspect, 
we adopt the following standard terminology: 

• an admissible closed arc (7, [a, 6]) in Vn+i is the restriction to a closed interval 
[a, 6] of an admissible section 7 : {c,d) — V„+i defined on some open interval 
(c,d) D [a,h]; 

• a piecewise differentiahle evolution of the system in the interval [tg, ^i] is a finite 
collection 

(7, [to,ti]) := {(7^''\ K-i,as]), s = l,...,N, to = ao < ai < ■ ■ ■ < aN = ti] 
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of admissible closed arcs satisfying the matching conditions 

7^ (as) = 7^'+'^ ys = l,...,N-l (2.38) 

Due to eq. (j2.38p . the image 7(i) is well defined and continuous for all to ^ 
t < ti, thus allowing to regard the map 7 : [toj^i] — *■ Vn+i as a section in a 
broad sense. The points 7(^0)) 7(^1) are called the endpoints of 7. The points 
Xs := 7(0^) , s = 1, . . . , — 1 are called the corners of 7. 

Consistently with the stated definitions, the lift of an admissible closed arc 
(7, [a, 6]) is the restriction to [a,b] of the lift 7 : {c,d) A, while the lift 7 of 
a piecewise differentiable evolution {(7^'''', [as_i,as])} is the family of lifts 7^'^^ 
each restricted to the interval [as-i,as]. 

With this definition, the image ■j{t) is well defined for all t 7^ ai, . . . ,aAr_i, 
thus allowing to regard 7 : [toj^i] — > ^ as a (generally discontinuous) section of 
the velocity space. In particular, due to the fact that the map i : A ^ ji(V„+i) 
is an imbedding of A into an affine bundle over V„+i , each difference 

[7],, = iil^'^'H^s)) -i{l^'Has)) ,s = l,...,N-l 

identifies a vertical vector in Tx^iVn+i) , henceforth called the jump of 7 at the 
corner Xg ■ 

In local coordinates, setting q^{'y^'^\t)) := g^*j,^(t), eqs. (|2.11|) . (j2.38|) provide 
the representation 

with [^^(7)] ^ := V'Ht^^^^H'^s)) ~ 4'^ i^^^Ho's)) denoting the jump of the function 
Vilit)) at t = a,. 

Going on with the generalization process, an admissible deformation of an ad- 
missible closed arc (7,(0,6]) is a 1-parameter family (75, [a(0) KOI) ! ICI < ^ 
of admissible closed arcs depending differentiably on ^ and satisfying the condi- 
tion (70, [a(0), 6(0)]) = (7, [a, 6]). Notice that the definition explicitly includes 
possible variations of the reference intervals [a(^),6(^)]. 

In a similar way, an admissible deformation of a piecewise differentiable evo- 
lution (7, [to,ti]) is a collection { (7^|\ [as_i(^), as(^)]) } of deformations of the 
various arcs, satisfying the matching conditions 

jfiasiO) = l^t'^iasiO) V|C| <e, s = l,...,iV-l (2.40) 

(s) 

Under the stated circumstances, the lifts 7^ and 7^ 1 respectively restricted 
to the intervals [a(^),6(^)] and [a5_i(^), as(^) ) are easily recognized to provide 
deformations for the lifts 7 : [a, 6] — >■ ^ and 7^*^ : [ag-i, as] — > A. 

Unless otherwise stated, we shall only consider deformations leaving the in- 
terval [toi^i] fixed, i.e. satisfying the conditions ao(C) = 1 o.n{0 = ti- No 
restriction will be posed on the functions as{S,) , s = 1, . . . , N — 1 . 
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Each curve Xs{S,) := 7^(as(^)) will be called the orbit of the corner Xs under 
the given deformation. 

In local coordinates, setting q^{'y^^\t)) = ip^^^{^,t) , the matching conditions 
read 

= <P(i+i)(e,«.(e)) (2.41) 

while the representation of the orbit Xs{S,) takes the form 

XsiO : t = asiO , = <^(!,)(C, (2-42) 

The previous arguments are naturally reflected into the definition of the in- 
finitesimal deformations. Thus, an admissible infinitesimal deformation of an ad- 
missible closed arc (7, [a, 6]) is a triple {a,X,/3), where X is the restriction to 
[a,b] of an admissible infinitesimal deformation of 7 : (c, d) — > Vn+i, while a, /3 
are the derivatives 



da 



(2.43) 

5=0 



expressing the speed of variation of the interval [a(^), b{^)] at ^ = 0. 

In a similar way, an admissible infinitesimal deformation of a piecewise differ- 
entiable evolution (7, [to,ii]) is a collection | - • • Os-i, X{s)i Os ■ ■ ■} of admissi- 
ble infinitesimal deformations of each single closed arc, with = ^■|g_Q, and, 

in particular, with oq = oat = whenever the interval [toi^i] is held fixed. 

Let us analyse the situation in detail. To start with, we notice that the quan- 
tities as, X^g-j are not independent: eqs. (j2.41|) imply in fact the identities 

9^(s) ^^das_ ^ ^V^Cs+i) ^^js+i) das 
dt di di dt di 

From these, evaluating everything at ^ = and recalling the relation between 
finite deformations and infinitesimal ones, we get the jump relations 



Moreover, the admissibility of each single infinitesimal deformation X^^s) es- 
quires the existence of a corresponding lift X(^s) = -'^(s) (^)^(.) + 
satisfying the variational equation (12.230 . 

Both aspects are conveniently accounted for assigning to each 7^'^) an (arbitrar- 
ily chosen) infinitesimal control : ^(7^'^^) — > ^(7^'^^) . In this way, proceeding 
as in § 2.3 and denoting by (■^)^(s) the absolute time derivative along 7^^^ induced 
by ^(s) we get the following 

Proposition 2.4 Every admissible infinitesimal deformation of an admissible evo- 
lution (7, [to,ii]) over a fixed interval [toj^i] is determined, up to initial data, by 
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a collection of vertical vector fields {^(s) = Y^-^ (af^)'y(s)}' s = 1, . . . , N and by 
N — 1 real numbers ai, . . . , a^-i through the covariant variational equations 



7 

completed with the jump conditions ^2.44^ - The lift of the deformation is described 
by the family of vector fields 

= +y(,), s=l,...,N (2.46) 

The proof is entirely straightforward, and is left to the reader. Introducing n piece- 
wise differentiable vector fields di, . . . ,dn along 7 according to the prescription 

= h(s) ( W-] V t G (as_i,a^) , s = l,...,N 

\OQ (t) 

eq. (I2.46P takes the explicit form 

^is) = hs) [H) j + y^s) = ^(1) + [-i^)^ (2.47) 

on each open arc 7^'^^ : {as-i,as) — > A. 

To discuss the implications of eq. (j2.45p we resume the notation ^^(7) for the 
totality of vertical vectors along 7. Notice that this makes perfectly good sense 
also at the corners 7(as)- 

We then define a transport law in ^^(7), henceforth called /i-transport, glu- 
ing /i(5)-transport along each arc (7^''^ [as_i,as]) and continuity at the corners, 
i.e. continuity of the components at t = . 

In view of Proposition 12.31 the /i-transported fields form an n-dimensional 
vector space Vh, isomorphic to each fibre 1^(7) |f • 

This provides a canonical identification of V{'~f) with the cartesian product 
[to,ii] X Vh, thus allowing to regard every section X : [to,ti] ^ ^(7) as a vector 
valued function X : [to,ti] ^ Vh- 

Exactly as in § 2.3, the situation is formalized referring V/j to a basis {e^^)} 

related to the basis (^^^ by the transformation 



7 

Given any admissible infinitesimal deformation {(X^^), [a^-i, as])}, we now 
glue all sections X^g-j : [as-i,as] ^(7*^*^) into a single, piecewise differentiable 
function X : [toj^i] ^ Vh, with jump discontinuities at t = expressed in 
components by eq. (I2.44p . For each s = 1, . . . , N this provides the representation 

/DX(g)\ dX"^ 
=X'^(t)e(,), [-^^ ^^=—^e^a) yte{as.,,as) (2.49) 
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In a similar way, we collect all fields Y(j,) into a single object Y , henceforth 
conventionally called a vertical vector field along 7. 

By abuse of language, we also denote by y = vector field along 

the open arcs of 7 defined by the prescription. 

y^(t) = as-i<t<as, s = l,...,N 

In this way, the covariant variational equation (|2.45p takes the form 



dt 



Y^e 



completed with the jump conditions 



1,...,A^-1 



(2.50) 



(2.51a) 



(2.51b) 



2.5 The variational setup 

A deeper insight into the algorithm discussed in § 2.4 is gained denoting by 53 
the infinite dimensional vector space formed by the totality of vertical vector fields 
Y = {Y(3) , s = 1, . . . , AT} along 7, and setting 2)3 := 2J M^-^ On account of 
eqs. ()2.51b .b). every admissible infinitesimal deformation of 7 is then determined, 
up to initial data, by an element (y, ai, . . . , aN-i) G 213. 

Let us now focus on the fact that, in the development of control theory, one 
is mainly interested in infinitesimal deformations A : [to, ti] ^ ^(7) vanishing at 
the endpoints. Setting X(to) = 0, eqs. ()2.51l a.b) provide the evaluation 

7 as<t 



X{t) = { I Y^e'^> ( ^ j dt-^as e^,'{as) ) e(,) (2.52) 




The vanishing of both A (to) and X{ti) is therefore expressed by the condition 

^^^^t^{^)dt - = (2.53) 

The left hand side of eq. ()2.53p defines a linear map T : 2B ^ V/i . According to 
the previous discussion, the kernel ker(T) C 2IT is isomorphic to the vector space 
of the admissible infinitesimal deformations vanishing at the end points of 7. 

Equally important is the nature of the inclusion T(213) C Vh- Depending 
on the latter, the evolutions of the system will be classified into normal, when 
T(2IJ) = Vfi and abnormal, when T(2B) C Vh- As we shall see, when applied to 
the extremals of an action functional, this terminology agrees with the current 
one (see, among others, \18\ |20] and references therein). The dimension of the 
annihilator (T(2]3))'^ C will be called the abnormality index of 7. 

In this connection, a useful characterization is provided by the following 
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Proposition 2.5 The annihilator (T(2IJ))'^ C Vj^ coincides with the totality of 
h-transported virtual 1-forms p = pid)^ satisfying the conditions 

P^(^) =0 ^ = l,...,r (2.54a) 



7 




Pi{as)[^p\J)]^^ = s = l,...,N-l (2.54b) 

Proof. In view of eq. (I2.53|) . the subspace (T(2n))'' C consists of the totahty 
of elements p = pa e^"^ = pa e^^^ Of' satisfying the relation 

V (y, ai, . . . , aAr_i) € 2n, clearly equivalent to eqs. (j2.54l a.b). □ 

On account of eqs. (j2.34l b). (j2.35l b). the condition of /i-transport of p along 
each arc 7^*^ is expressed in coordinates as 

the cancellation arising from the requirement ()2.54l al. 

The content of Proposition 12.51 is therefore independent of the choice of the 
infinitesimal controls : ^(7*^*^) ^(7^'^)) . 

Remark 2.2 According to Proposition \2. 51 the abnormality index of a piecewise 
differ entiahle section 7 cannot exceed the abnormality index of each single arc 
7^^) . Thus, for example, if one of the arcs is normal, 7 is necessarily normal. 
More generally, due to the additional restrictions posed by eqs. \2.54-t >) and by the 
continuity requirements [p]as = 0, an evolution may happen to be normal even if 
all its arcs 7*-^^ are abnormal. Typical examples are: 

• Vn+i = M X £"2, referred to coordinates t,x,y. Constraint: iP' + ip' = f^. 
Imbedding A — > ji(V„+i) expressed in coordinates as x = vcosz, y = vsinz. 
Piecewise differentiable evolution 7 consisting of two arcs: 

7^^) : X = 0, y = vt to<t<0 
^(2) : x = vt, y = < t < ti 

Eq. Ii2.54[ i) admits h-transported solutions p^^^ = aCo^ along 'j^'^^ and p^^ = /?cj^ 
along 7(2), V a,/3 G M. Both arcs are therefore abnormal. Notwithstanding, 7 is 
normal, since no pair p^^\ p^^ matches into a continuous non-null virtual 1-form 
along 7. 

• V„+i=IRxi?2- Coordinates t, x,y . Constraint: v'^ x = {ij^ — t^)"^ . Imbedding 
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^ ~^ Ji(Vn+i) expressed in coordinates as x = v ^(z^ — a^t^)^, y = z. Piecewise 
dijferentiable evolution 7 consisting of two arcs: 

7(^)1 x = 0, y = ^a{t'^ -t*^) to<t<t* 

7^^^: x = -^{t^ -t*^), y = t* < t < ti 

{t* 7^ j. -Eg. ^2. 54P -) admits h-transported solutions of the form p = atv^ along 
the whole of 7 . Both arcs 7^"^^ , 7^^^^ are therefore abnormal. Notwithstanding, 7 
is normal, since no solution satisfies condition {2. 54p )- 

As a concluding remark, let us finally observe that, although geometrically 
significant, the arguments discussed so far provide only a partial picture of the 
situation. In a variational context, in fact, what really matters is not the totality of 
admissible infinitesimal deformations vanishing at the end points — here identified 
with the kernel of the map T : 213 — > V/j — but the (possibly smaller) subfamily X 
of infinitesimal deformations tangent to admissible finite deformations with fixed 
end points. 

The linear span of X, henceforth denoted by A(7) , will be called the variational 
space of 7. The evolutions of the system will be classified into ordinary, when 
A(7) = ker(T) and exceptional, when A(7) C ker(T). A hierarchy between the 
various typologies is provided by the following 

Proposition 2.6 The normal evolutions form a subset of the ordinary ones. 
The result is proved in Appendix B. In this connection, see also |24]. 

3 Calculus of variations 
3.1 Extremals 

Let us now come to the central problem of control theory. Let L G 5^(^) denote a 
difFerentiable function on the velocity space A, henceforth called the Lagrangian. 

Also, let (7, [to, ^i]) (7 for short) denote an admissible piecewise difFerentiable 
evolution of the system, defined on a closed interval [toi^i] C M. 

Indicating by 7 the lift of 7 , define the action functional 

/[7]:= Ldt:=Y. / {i^'^)\L) dt (3.1) 

•'l s = l O.S-1 

Definition 3.1 An admissible evolution 7 is called an extremal for the functional 
113. 1\) if and only if, for all admissible deformations 7^ = {(7^|'', [0^5-1(0) ) } 

with fixed endpoints, the function 
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\L)dt 



is stationary at ^ = . 

The content of Definition 13.11 is formalized denoting by X(^g) the infinitesimal 
deformation associated with each single 7 c*'' . Recalling eq. (12.470 as well as the 



definition a, = ^ 



1^=0 



this provides the evaluation 



dl 

d^ 



N 

= E 

«=0 s=l 



N 

E 

s=l 



X(^s){L)dt + \asL{^'^'Has))-as-iL{j^'\as-i)) 



a.s-1 



(3.2a) 



On account of the assumption oq = oat = , denoting by 

the jump of the function L{'j{t)) at t = ag, eq. ()3.2l a) may be concisely written as 



dl_ 



/as , r)T \ 

[xlg{d.{L)+Y^%-^)dt - Y^^^a. (3-2b) 



Eq. (j3.2l b) is further elaborated introducing N virtual 1-form A^'') = {t) 
(one for each arc 7'-^'' ) satisfying the transport law 



Dt 



diL 



as well as the matching conditions 



A(^+i)| 



v{s) 



s = l, 



(3.3a) 



A^- 1 



(3.3b) 



Once again, for notational convenience, we collect all A^'^^ into a continuous, 
piecewise differentiable section A : [^Oi^i] — ^ ^*(7) according to the prescription 



A(t) = A(^)(t) VtG[a,_i,a,] 



(3.4) 



On account of eqs. (|3.3b .b). A is then uniquely determined by L, up to initial 
data at t = t^. 
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Taking the covariant variational equation (j2.45p as well as the duality relations 



y(s) 



6^ into account, by eq. ()3.3l a) we get the expression 



whence also 



Dt 



IE 



Dt 



P 



is) 



a.s-1 



^(1) Ul) dt 



X 



is)P 



as-1 



Summing over s, restoring the notations ()2.50p . (j3.4p and recalling eqs. (j2.44p . 
()3.3b ) as well as the conditions X(to) = -^(^i) = 0, this implies the relation 



N 

E 

s=l 



X(s) ^^{L) dt 



to 



N-l 



L J a. 



s=l 



In this way, omitting all unnecessary subscripts, eq. ()3.2l b) gets the final form 



dl 



to 



dL 



df_\ 



N-l 



Y^dt + Y,as[p;{t)iP'{^)-m 



s=l 



(3.5) 



In the algebraic environment introduced in § 2.5, the previous discussion is 
naturally formalized regarding the right hand side of eq. ()3.5p as a linear functional 
diry : 2ir ^ M on the vector space 2IJ = 53 © M^~^. A necessary and sufficient 
condition for 7 to be an extremal for the functional (13. ip is then the vanishing 
of dl^ on the subset X C 2)3 formed by the totality of elements Y,ai,... ,aj\f-i 
arising from finite deformations with fixed end points. By linearity, the previous 
condition is mathematically equivalent to the requirement 



A (7) C ker(d/^ 



(3.6) 



with A(7) = Span(X) C ker(T) denoting the variational space of 7. 

As we shall see, eq. (13. 6p provides an algorithm for the determination of all the 
extremals of the functional (j3.ip within the class of ordinary evolutions. 

The exceptional case is considerably more complicated, due to the lack of an 
explicit characterization of the space A (7) in terms of the local properties of 
the section 7. In this respect, the simplest procedure and, quite often, the only 
available one, is checking eq. (13. 6p separately on each exceptional evolution. 

In what follows we shall adopt an intermediate strategy, namely, rather than 
dealing with eq. (13. 6p we shall discuss the implications of the stronger requirement 



ker(T) C ker(d/^) 



(3.7a) 
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According to the classification introduced in § 2.5, the latter is necessary and 
sufficient for an ordinary evolution 7 to be an extremal of the functional (j3.ip . 
but merely sufficient for an exceptional evolution to be an extremal. 

By elementary algebra, the requirement (j3.7l a) is equivalent to the existence of 
a (possibly non-unique) linear functional i^T : V/j — > M satisfying the relation 



T 

2n Vh 




(3.7b) 



Setting K = Kae^"'\ and recalling eqs. p.53p . (13. 5p . the requirement p. 71 b) is 
expressed in components as 



to \dz^ 



dtp' 



N-l 



Y^dt + Y.a,[p,{tW{^)-m 

N-l 



dz 



V^(7) 



By the arbitrariness of y, ai, . . . , a^-\ , the latter condition splits into the system 
dL 







p,+Kae^f)ij\l)-m) 







^ = l,...,r 
s = l,...,iV- 1 



(3.8a) 
(3.8b) 



Collecting all results, and recalling Propositions 12.5 1 12.61 we conclude 



Theorem 3.1 Given an admissible evolution 7, let p(7) denote the totality of 
piecewise differentiable virtual Informs A = pi{t)uj'^ along 7 satisfying eqs. h3.3\a .h). 
{3.4^ as well as the finite relations 



dz^ 



dL 

dz^ 



and the matching conditions 



A = l,...,r (3.9a) 



s = l,...,N -1 (3.9b) 



Then: 

a) the condition p{'y) ^ 'i) is sufficient for 7 to he an extremal for the functional 

(EIP; 
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b) if ^ is an ordinary evolution, the same condition is also necessary for 7 to 
be an extremal; 

c) J is a normal extremal, namely an extremal belonging to the class of normal 
evolutions, if and only if the set p{'y) consists of a single element. 

Proof. In view of eqs. (|3.5p . (|3.9l a.b). the ansatz A G p(7), if allowed, implies 
= for all admissible infinitesimal deformations vanishing at the end points 
of 7. Assertion a) is then a direct consequence of Definition 13. 1[ 

In particular, according to our previous discussion, if 7 is an ordinary extremal, 
for any continuous virtual 1-form A = piOJ^ obeying the transport law (jS.SI a) there 
exists at least one /i~transported 1-form K = KaC^"^ satisfying eq. (I3.8l a.b). The 
sum \ + K = (^pi +Ka e^f^ ) a)* is then automatically in the class ^(7), thus proving 
assertion b). 

Finally, as pointed out in § 2.4, the normal evolutions form a subclass of 
the ordinary ones, uniquely characterized by the requirement (T(2IJ))'^ = {0}. 
Therefore, according to assertion b), a normal evolution 7 is an extremal if and 
only if the class p(7) is nonempty. Moreover, by eqs. (|3.3l a). (|3.8l a). if A, A' is 
any pair of elements in the class p{"f) , the difference A — A' is automatically an 
/i-transported 1-form satisfying eqs. (|2.54l a.b). By Proposition 12.51 this implies 
A - A' G (T(2n))° ^ A = A', thus establishing assertion c). □ 

In view of eqs. (j2.34l bl. (j2.35l b). for any A G p(7) the transport law (j3.3l a) 
simplifies to 




the cancellation being due to eq. ()3.9l a). Exactly as it happened with Proposi- 
tion 12.51 all assertions of Theorem 13.11 have therefore an intrinsic meaning, irre- 
spective of the choice of the infinitesimal controls : V{j^'^^) ^(7^'^)) . 

The previous arguments provide an algorithm for the determination of the 
ordinary extremals of the functional ()3.ip . relying on 2n + r equations 



dq' 
'dt 



ip\t,q\z^) (3.10a) 



dpi di)^ dL /QinKA 

-dt^W'^^W ^'-'"^^ 

d^p' dL 

P^^ = ^ (3-lOc) 

for the unknowns q^{t) ,pi{t) , z^(t)., completed with the continuity requirements 

W]a. = = [P^^'-L]^^ =0 . = l,...,iV_i (3.11) 

As already pointed out, all equations are independent of the choice of the 
infinitesimal controls, and involve only the "true" data of the problem, namely the 
Lagrangian L and the constraint equations (I2.10p . In particular: 
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• the algorithm (j3.10l a.b.c). (jS.lip is invariant under arbitrary transformations 
of the form 

f{t, q^, . . . , q"-) being any differentiable function over Vn+i ; 

• the last pair of equations (j3.1ip extend to the ordinary evolutions the well 
known Erdmann- Weierstrass corner conditions of holonomic variational cal- 
culus [g do] . 

3.2 Lagrange multipliers 

For completeness, we discuss the relation between Theorem 13.11 and the traditional 
approach to constrained variational calculus based on the so called Lagrange mul- 
tipliers method. To this end, we consider a situation in which: 

• the behaviour of the system is determined by a "free" lagrangian C, viewed 
as a function on ji(V„+i), with local expression C = C{t,q^,q^); 

• the presence of the constraints restricts the mobility of the system to a sub- 
manifold A ji(Vn+i), implicitly represented as 

=0 a = l,...,n-r 

with g"" G 5(ii(V„+i)) and rank H^^^r^r^ \\=n-r. 

The geometric setup developed in § 3.1 is recovered regarding the intrinsic 
lagrangian as the pull-back of the extrinsic one. In coordinates, the resulting 
state of affairs is summarized into the equations 

L(t,g\z^) = Cit,q\^l^\t,q\z^)) (3.13a) 
g^t,q\^l^\t,q\z^)) = (3.13b) 

These, in turn, imply the relations 

dL _ dC dC dijj^ dL _ dC dt/j^ 

dq^ dq^ dq^ dq^ ' dz^ dq^ dz^ 

dg^ dg^ dip^ _ ^ dg^dip^ _ ^ 

dq^ dq^ dq^ ' dq^ dz^ 

Let us now consider the product manifold £ := Vn+i x M"^^, referred to fiber 
coordinates t,g*,Acr. The natural projection £ — Vn+i makes £ into a vector 
bundle over Vn+i and therefore also into a fiber bundle over M. We regard £ as 
the event space of a fictitious unconstrained system, glue the functions C, g'^ € 
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into a single function £ := £ + A^i?'^ € and adopt £ as a 

(singular) lagrangian for the newer system. 

In this way, proceeding as in § 3.1, to every continuous, piecewise differentiable 
section /i : [toi^i] — we associate the action integral 



^ '-as , 



.8 = 1 



I[fi]:= jCdt:=Y, C{t,q\q') + \„g%t,q\ci) dt (3.15) 



fi denoting the lift of ;U to a section fi : [toi^i] ^ 3i{£)- We can then state 

Theorem 3.2 The projection tti : £ ^ Vn+i sets up a 1-1 correspondence be- 
tween extremals of the functional hS.l^^) and extremals of the functional hS. 1\) sat- 
isfying condition a) of Theorem \3.1\ 



Proof. According to elementary (holonomic) variational calculus, the extremals 
of the functional ()3.15p are determined by the Euler-Lagrange equations 

g%t,q\q') = Q (3.16a) 

dt'd^ ~ dq^ ~ y ■ ) 

completed by the Erdmann-Weierstrass corner conditions, asserting the continuity 

of the "momenta" ^ and of the "hamiltonian" {q^ ~ ^) along fi. 

On the other hand, eqs. ()3.14l a.b) and the definition of C imply the relations 

dq'^ dq" dq^ dq^ ' dz^ dq^ dz^ 

From these, taking eqs. (I3.13l a.b). (I3.14l a.b) into account, we conclude: 

• if /i : = q^{t), = \a{t) is an extremal for the functional ()3.15p . the 
functions q\t) and pi{t) := /i*(^) satisfy eqs. (fXTUl a.b.c). (fTTTT) . There- 
fore, the class p('/ri • /i) associated with the section vri • ^ : M ^ Vn+i is 

dC 



non-empty, since it contains at least the 1-form A = /x*(^) cu 



conversely, let q^ = q^{t) , pi = Pi{t) be any solution of eqs. ()3.10l a.b.c). ()3.1ip : 
then eqs. (j3.1Ul a). (j3.13l b) automatically imply eq. (j3.16l a) . Moreover, on 
account of the identification (13. 131 a). eq. (I3.10fc ) reads 

dijj' dL ( dC\ dip' 



dz^ dz^ \ * dq^ J dz^ 



ensures that the linear system 



In view of the assumption rank || g(gi!.^gn) || = n — r, the second equation ()3.14b ) 



dg"^ dC 
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can be uniquely solved for the functions = Ao-(t). 
Together with eqs. (j3.14l a.b). this implies the further relation 



d^lj^ dL _ / dC ^ dg''\dip'' dL _ dC ^ dg" 
dq^ dq^ \ dq^ ^ dq'^ J dq^ dq^ dq^ ^ dq^ 

Eqs. p.lSl a.b) allow to cast eq. ()3.10l b) in the equivalent form 

dt dq^ dq^ dt \ dcf " dq^ J dq^ " dq^ 

reproducing the content of eqs. ()3.16l b). with C = C + X^g" ■ 
Finally, on account of eq. (j3.18l a). eqs. (jS.lip are easily seen to coincide with the 
Erdmann-Weierstrass conditions along the curve = q^{t), A^- = Ao-(t). Every 
solution = q^{t),pi = Pi{t) of the system ()3.10l a.b.cl. (jS.lip is therefore related 
to a unique extremal of the functional ()3.15p by the projection algorithm. □ 



3.3 Pontryagin's "maximum principle" revisited 

From a mathematical viewpoint, eqs. (|3.1Ul a.b.c). form a set of 2n + r equations 
for the 2n + r unknowns q''{t), z'^{t) , pi{t) . In this respect, they have a natural 
setting in the geometrical environment provided by the contact bundle C{A) . 

As pointed out in Appendix A, the latter is a vector bundle over A, identical 
to the pull-back of the phase space l^*(V„+i) through the commutative diagram 

C{A) V*{Vn+l) 

C 

A Vn+l 

Accordingly, we refer C{A) to fibered coordinates t,q'',z^,pi , with t,q^,z'^ coor- 
dinates in A and t,q^,pi coordinates in V*{Vn+i) ■ 

The advantage of the environment C{A) comes from the presence of a distin- 
guished Liouville 1-form Q, expressed in coordinates as 

e = piCb' = pi{dq' - ip'dt) (3.19) 

(For simplicity, we are not distinguishing between covariant objects in A and their 
pull-back in C{A), namely we are writing ■0* for C*('0*), uj^ for (^*(uj^) etc.) 

By means of Q, every Lagrangian L S d{A) may be lifted to a 1-form -Sl over 
C{A) according to the prescription 

T^L ■■= Ldt + e = {L-pitp')dt+pidq' := -Hdt+pidq' (3.20) 

The difference TC := piip^ — L is known in the literature as the Pontryagin 
Hamiltonian. Needless to say, 7i is not a Hamiltonian in the traditional sense, 
but a function on the contact bundle. 
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To understand the role of the -Ol we focus on the fibration C{A) V„+i given 

by the composite map v := tt ■ k. A piecewise differentiate section (cr, [to,ii]) 
consisting of a finite family of closed arcs 
As) . 



1, 



, A'', to = aQ < ai < ■ ■ ■ < aN = ti 



will be called v-continuous if and only if the composite map f • o" is continuous, 
i.e. if and only if a projects onto a continuous, piecewise differentiable section 

v-a : [to)*i] Vn+i- A deformation = {{(^^f , [o-s-iiO^'^siO]) } ^i^l similarly 
be called u-continuous if and only if all sections are ^-continuous. A necessary 
and sufficient condition for this to happen is the validity of the matching conditions 
()2.40p . synthetically written as 

s = l,...,A^-l (3.21) 



lim V ■ cr^{t) 



lim V ■ (t) 



A f-continuous deformation is said to preserve the end points of f • o" if and 
only if f • fj^ is a deformation with fixed endpoints. A vector field along a tangent 
to the orbits of a u-continuous deformation is called an infinitesimal deformation. 

Notice that, since the stated definitions do not include any admissibility re- 
quirement for the sections v ■ a^, the only condition needed in order for a vector 
field X* (^)o. + r'^(g^)^ + nj(^)^ to represent an infinitesimal deformation of 
a is the consistency with the matching conditions (j3.2ip . expressed in components 
by the jump relations 



lim 



X' + a. 



dq^ 
It 



lim 



X' + a, 



dq^ 



1,...,N-1 (3.22) 



' das \ 



with Os 

By means 'di we now define an action integral over C{A) , assigning to each 
u-continuous section a : q^ = q^ (t) , = (t) , Pi = Pi (t) the real number 



to 



dt 



(3.23) 



For any u-continuous deformations preserving the end points of u • a we have 
then the relation 



dl[(7^] 

N 
s=l 



lim 



f'^ \fdl_dn\ _(dm.dn\ 

k \\dt \dt ^ dqO 



dn^A 

dz^ 



dt + 



dq' 
'dt 



n] +PrX' 



lim 



dt 



n] +p^x' 



From the latter, taking eqs. (|3.22|) and the conditions X*(to) = X^{ti) = 
into account, we conclude that the vanishing of ^|^_q under arbitrary defor- 
mations of the given class is mathematically equivalent to the system 



d(l_ 
'dt 



dn 

dpi 



dpi 
dt 



dq^ ' 



dz^ 



(3.24a) 
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completed with the continuity conditions 



[P^]as=Mas=^ . = l,...,iV-l (3.24b) 

where, as usual, we are denoting by [f]as the jump of the function f{t) at t = ag. 

In view of the definition of the Pontryagin Hamiltonian Tl, eqs. (|3.24l a.b) are 
easily seen to reproduce the content of eqs. ()3.10l a.b.c). (13. lip , the continuity 
requirement [q'^]as = being implicit in the definition of a. 

As far as the ordinary extremals are concerned, the original (constrained) vari- 
ational problem in the event space is therefore equivalent to a free variational 
problem in the contact manifold. This is precisely the essence of Pontryagin's 
maximum principle. 

As a further comment on eqs. (13.240 . let us digress on the special situation 
determined by the ansatz L = 0. Under the stated circumstance, the functional 

Xo[a] ■■= = P^ {% - V-^) dt (3.25) 

is an intrinsic attribute of the manifold C{A) , entirely determined by the Liouville 
1-form ()3.19p . Its role is clarified by the following 

Proposition 3.1 Let 7 : [toj^i] Vn+i denote any continuous, piecewise differ- 
entiable section. Then: 

a) ^ is admissible if and only if the functional 113. 25\) admits at least one ex- 
tremal a projecting onto 7, i.e. satisfying v ■ a = j ; 

h) for any such 7, the totality of extremals of Iq projecting onto 7 form a finite 
dimensional vector space over M, with dimension equal to the abnormality 
index of 7 . 

Proof. For L = 0, eqs. (j3.24l a.b) reduce to 

(3.26a) 
(3.26b) 
(3.26c) 

s = l,...,iV-l (3.26d) 

Eq. (13. 261 a) is the admissibility requirement for the section v ■ a . Due to this 
fact, if an extremal a of the functional (j3.25p satisfies f • cr = 7, its projection C, ■ o 
under the map C '■ C(^) — > A coincides with the lift 7 : [toi^i] A. 

For any admissible 7, the extremals projecting onto 7 are therefore in 1-1 
correspondence with the solutions Pi(i) of the homogeneous system (I3.26l b.c.d). 



dpi di)'^ 



p 
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with the functions z^{t) regarded as given. On the other hand, according 

to Proposition 12.51 . eqs. (|3.26l b.c.d) are precisely the relations characterizing the 
totality of virtual 1-forms pi{t)id^ belonging to the annihilator (T(2IT))'^. 
Both assertions a) and b) follow easily from this fact. □ 

In the language of § 2.5, Proposition 13.11 asserts that a section 7 : [to,ii] Vn+i 
describes a normal evolution of the system if and only if the functional ()3.25p 
admits exactly one extremal projecting onto 7, namely the one corresponding to 
the trivial solution pi{t) = 0. If the extremals projecting onto 7 are more than one, 
7 represents an abnormal evolution; if no such extremal exists, 7 is not admissible. 
Returning to the action integral (j3.23p we can now state 

Proposition 3.2 The totality of extremals of the functional ^3. 23\) projecting onto 
a section 7 : [tQ, ti] V„+i is an affine space, modelled on the vector space formed 
by the extremals of the functional 13. 25\) projecting onto 7. 

The proof, entirely straightforward, is left to the reader. 

The previous arguments provide a restatement of Theorem 13.11 in the environ- 
ment C{A) . In particular, it is worth remarking that, in general, the projection 
algorithm a ^ v ■ cr , applied to the totality of extremals of the functional ()3.23p . 
does not yield back all the extremals of the functional (j3.ip . but only a subclass, 
wide enough to include the ordinary ones. The missing extremals may be ob- 
tained determining the abnormal evolutions by means of Proposition 13. H finding 
out which ones have an exceptional character, and analysing each of them individ- 
ually. 



3.4 Hamiltonian formulation 

As pointed out in § 3.2, all ordinary extremals of the functional (j3.ip are projections 
of extremals of the functional (I3.23p . Let us analyse the implications of this fact. 

To this end, temporarily leaving aside all aspects related to the presence of 
corners, we observe that a differentiable curve a in C{A) is at the same time a 
section with respect to the fibration C{A) M and an extremal for the functional 
(|3.23p if and only if its tangent vector field Z := satisfies the properties 

{Z,dt) = l, ZJd'dL = (3.27) 

On account of eq. (j3.20p . at any G C{A) a necessary and sufficient condition 
for the existence of at least one vector Z € T^(C(^)) satisfying eqs. ()3.27p is the 
validity of the relations 

1^)^ = (3-28a) 

Points ? at which eqs. (j3.27p admit a unique solution Z will be called regular 
points for the functional (I3.23p . In coordinates, the regularity requirement is 
expressed by the condition 
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In view of eq. (j3.28l b). in a neighborhood of each regular point eqs. (|3.28h ) may 
be solved for the z'^'s, giving rise to a representation of the form 

= z^{t,q\...,q^,pi,...,pn) (3.29) 

The regular points form therefore a (2n + l)-dimensional submanifold S C{A), 
locally diffeomorphic to the phase space V*{Vn+i)- 

Inserting eqs. (j3.29p into the first pair of relations ()3.24l a) gives rise to a system 
of ordinary differential equations in normal form for the unknowns q^{t),Pi{t) . The 
algorithm is readily implemented denoting by H := the pull-back of the 

Pontryagin Hamiltonian to the submanifold S, expressed in coordinates as 

H = n{t,q\z^{t,q\p,),pr) = Pki^Ht,q\z^{t,q\Pi)) - L{t,q\ z^{t,q\pi)) 
In view of eqs. ()3.20p . ()3.28l a) we have then the identifications 

= = (3.30a) 

dpi dpi ydz^ dpi 

m _ dn dp^dz^ _ dL_ 
dq^ dq^ /6z^ dq^ dq^ dq^ 

On account of these, the first pair of equations (13. 241 a) takes the final form 
dq' dH 

^=-^ (3.31b) 
dt dq' ^ ' 

The original constrained Lagrangian variational problem has thus been reduced 
to a free Hamiltonian problem on the submanifold j : 5 — > C(^), with Hamilto- 
nian R{t^ g^, . . . , (f' ^p\ , . . . ,Pn) identical to the pull-back H = j*{'H). Conversely, 
setting H = j*(7i), the inverse Legendre transformation = together with 
eq. ()3.30l a). yields back the constraint equations g* = ip''{t,q^, z^) . Once again, all 
this is in full agreement with Pontryagin's principle. 

A f-continuous extremal of the functional (j3.23p consisting of a finite family 
of closed arcs a^'^^ : [0^-1, «s] C{A), each contained in (a connected component 
of) the submanifold S will be called a regular extremal. 

Singular extremals, partly, or even totally lying outside S may also exist. In 
fact, while eq. ()3.28l a) is part of the system ()3.24l a.b). and must therefore be 
satisfied by any extremal, the requirement (j3.28l b) has only to do with the well- 
posedness of the Cauchy problem for the subsystem ()3.24l a). 

On the other hand, by construction, the Hamilton equations (j3.31l a.bl deter- 
mines only the regular extremals. The singular ones, when present, have therefore 
to be dealt with directly, looking for solutions of eqs. (j3.24l a.b) not arising from 
a well posed Cauchy problem. In principle, this could be done extending to the 
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non-holonomic context the concepts and methods commonly adopted in the study 
of singular Lagrangians [57]. The argument is beyond the purposes of the present 
work, and will not be pursued. 

To complete our analysis, let us finally discuss the role of eqs. (|3.24l b) in the 
study of corners. To this end, we consider the cotangent space T*{Vn+i), referred 
to local fibered coordinates t,q^,po,Pi, and denote by 

^ := PQdt + pidq' (3.32) 

the corresponding Liouville 1-form. We also recall (see Appendix A) that the 
equivalence relation (|A.6p determines afibration T*{Vn+i) — > ^*(Vn+i), expressed 
in coordinates as {t,q^,pQ,pi) {t,q^,pi). 

Making use of the 1-forms (j3.20p . (|3.32|) let us now introduce a morphism 

CiA) ^T*{Vn+i) fibered over V*{Vn+i), based on the prescription 

In coordinates, we have the explicit representation 

t = t, q' = q\ p^ = Pi , po = -n{t,q\pi, z^) (3.33) 

The content of eqs. (j3.24l b) is then summarized into the following 

Proposition 3.3 For any v-continuous extremal (o", [to,ii]) of the functional 
113. 23\) . the composite map ■ a : [tQ,ti] — > T*(Vn+i) is necessarily continuous. 

The previous arguments provide a simple characterization of the jumps that 
may possibly occur along a regular extremal a : {[a^'^\ [as^i,as])} . To this end 
we observe that the restriction of the map (|3.33|) to the submanifold S C C{A) 
determines an immersion : S ^ T*(Vn+i) and that, as already pointed out, at 
each ? e 5 there exists, locally, one and only one differentiable extremal of the 
functional (I3.23P through 

On the other hand, by Proposition 13.31 for each s = 1, . . . , — 1 , the arcs a^^^ 
and (t(*+^) are related by the condition ^((j*^*)(as)) = ^ [a^'^~^^\as)) . From this 
it is readily seen that the admissible discontinuities of a or, what is the same, the 
admissible corners in the projection j := v ■ a : [to,ti] Vn+i may only occur at 
points in which the immersion ^' : 5 — > T*(V„-|-i) is not injective. 

A Geometry of the velocity space 

(i) Given the event space V„+i , let ji(Vn+i) and V{Vn+i) respectively denote the 
first-jet bundle and the vertical bundle associated with the fibration t : Vn+i — > K. 
Both spaces ji(Vn+i) and V{Vn+i) may be viewed as submanifolds of the tangent 
bundle T(Vn+i) according to the identifications 

ji(V„+i) = {z\ze r(V„+i) , {z,dt) = 1} (A.la) 

V{Vn+i) = {v \ veT{Vn+i), {v,dt)=0} (A. lb) 
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Eqs. (jA.ll a.b) point out the nature of ji(Vr!,+i) as an affine bundle over Vn+i, 
modelled on the vertical space V{Vn+i) [221 [Ml El IS] • 

Given any local fibered coordinate system t, over V^+i , we respectively refer 
ji(Vn+i) to jet-coordinates t,q^,q^ and V{Vn+i) to coordinates t,q^^u^. In terms 
of these, the content of eqs. ()A.ll a.b) is summarized into the relations 

- + q\z)--] V^Gii(V„+i) (A.2a) 

^ = ^'(^) ^ ^ e ViVn+i) (A.2b) 

Eq. ()A.2l a) allows to set up a vector-bundle homomorphism 

T{h{Vn+l)) V{Vn+l) 



ji(V„+i) — ^ V„+i 

assigning to each X G r2(ji(Vn+i)) a vertical vector ^{X) € 1/^(2) C^'n+i) according 
to the prescription 

v{X) := 7r,(X) - z (A.3) 

In coordinates, introducing the notation 

J := dq' - q'dt (A.4) 

eq. (|A.3p takes the explicit form 

(ii) The dual of the vertical bundle, henceforth denoted by V*{Vn+i) Vn+i, 
is called the phase space. In view of eq. (lA.lb ). the latter is canonically isomorphic 
to the quotient of the cotangent space T*(Vn+i) — ^ V„+i by the equivalence 
relation 

o" ~ cr -4=^ < (A.6) 

For simplicity, we preserve the notation ( , ) for the pairing between ViVn+i) 
andl/*(V„+i). 

Every local coordinate system t,q^ in Vn+i induces fibered coordinates t,q^,pi 
in y*(V„+i) , with pi{a) :=(&, (^) ) V a E V*{Vn+i) ■ 



7r(<3-) 



(iii) Let V{ji{Vn+i)) ^ ji{Vn+i) denote the vertical bundle over ji(Vn+i) relative 
to the fibration vr : ji(V„+i) V„+i . Given any jet coordinate system t, q*, in 
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ii(Vn+i), we refer V{ji{Vn+i)) to fibered coordinates t,q'^,q'^,u^ , with uniquely 
defined by the requirement 

The affine character of the fibration ji(V„+i) ^ Vn+i provides a canonical 
identification of F(ji(V„+i)) with the pull-back of the vertical bundle V{Vn+i) ■ 
This gives rise to a vector bundle homomorphism 

V{jl{Vn+l)) V{Vn+l) 

C TT (A.7a) 

ji(V„+i) — ^ V„+i 

expressed in coordinates as 

The result is well known from non-holonomic mechanics (see e.g jl6 t ll7 t [22 | l26j 
and references therein). A sketchy proof may be traced as follows: for each z G 
ji(Vn+i) , the fiber := 7r~^(7r(z)) through z is an affine submanifold of ji(Vn+i) , 
modelled on the vertical space K-(z)C^'n+i) • Every pair {z,v) , v G K-(2)C^'n+i) is 
therefore an "applied vector" at z in S^, i.e. an element of the tangent space 
Tz{T,z) ■ On the other hand, by definition, Tz{Tiz) is isomorphic to the vertical 
space VziJiiVn+i)) ■ By varying z we conclude that the totality of pairs {z,v) G 
ii(Vn+i) xv„+i V{Vn+i) is in 1-1 correspondence with the points of l^(ji(Vn+i)) , 
thus establishing the result. 

(iv) The pull-back of the space V*{Vn+i) through the map ji(Vn+i) — Vn+i 
determines another important space C(ji(Vn+i)), henceforth referred to as the 
contact bundle. Once again, the situation is conveniently summarized into a com- 
mutative diagram 

C(jl(V„+l)) V*{Vn+l) 

c| |. (A.8) 

Jl(V„+l) ^ Vn+1 

Notice that, by construction, C {ji{Vn+i)) is at the same time a vector bundle 
over jiiVn+i) and an affine bundle over V*{Vn+i)- 

The manifold C(ji(Vn+i)) will be referred to coordinates t,q'^,q^,pi related in 
an obvious way to the coordinates t,q^,q^ in ii(Vn+i) and t,q^,pi in V*{Vn+i) ■ 
Every a G C(ji(V„+i)) will be called a contact 1-form over ji(Vn+i)- 

In view of the stated definition, a contact 1-form a is essentially a pair (z, a) G 
ji(V„+i) X T/*(V„+i), with a G V*..{Vn+i)- Now, by eq. i^EM, every such pair 



34 



determines a linear functional on the tangent space Tz{ji{Vn+i)) according to the 
prescription 

{a,X) := {a,u{X)) V X G r,(ii(V„+i)) (A.9) 

In coordinates, recalling eq. ()A.5p . the definition of pj(o') and the identification 
Pi{cr) = Pi{a), eq. (jA.Qh takes the explicit form 

Prom the latter it is easily seen that the knowledge of the functional ()A.9p is 
mathematically equivalent to the knowledge of a . 

The contact bundle C(ji(Vn+i)) is therefore identical to the vector subbundle 
of the cotangent space T*(ji(Vn+i)) locally generated by the 1-forms ()A.4p . while 
the coordinates pi coincide with the components involved in the representation 



p,(a)c^*l^(,) VaeC(ii(V„+i)) (A.IO) 



An important attribute of the contact bundle is the presence of a distinguished 
Liouville 1-form 0, assigning to each a G C{ji{Vn+i)) the linear functional G 
T*[C{ji{Vn+i))] defined by the prescription 

{Qla,Y) := {a,UY)) V Y G r4C(ji(V„+i))] (A.ll) 

In coordinates, using the same notation for the 1-forms ()A.4p and for their 
pull-back on C{ji{Vn+i)) , eqs. (jA.lOp . (jA.lip yield the representation 

e = piuj' = pi{dq' - q'dt) (A.12) 

(v) In the presence of non-holonomic constraints, let A ii(Vn+i) denote the 
submanifold of ji(Vn+i) formed by the totality of admissible velocities. 

As in § 2.1, we refer A to fibered coordinates t, q^, z^, and represent the imbed- 
ding i : A^ ji{Vn+i) locally as 

q' = i^\t, g\ . . . , g", z\ . . . , z'') i = l,...,n 

The concepts of vertical vector and contact 1-form are immediately extended 
to the submanifold A : as usual, the vertical bundle V{A) is the kernel of the 
push-forward T{A) T(Vn+i), while the contact bundle C{A) is the pull-back 
on A of the contact bundle C(ji(Vn+i)), as expressed by the commutative diagram 

C{A) C(ji(V„+i)) 



c 



C (A.13) 



A Jl(Vn 



+1; 
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The manifolds V{A) , C{A) will be respectively referred to coordinates t, q\ z^, 
and t, q\ z^,pi . In this way, setting 



(A.14) 



we have the representations X = ^^{X) (^-^x) i^(^x)^ f'" = Pi(^) '^*|C(ct) ^X (^V{A), 
aeC{A). 

The pull-back := i*{Q) of the Liouville 1-form of C(ji(Vn+i)) will be called 
the Liouville 1-form of ^. In coordinates, eqs. ([02]) . (TOUl provide the repre- 
sentation 

Q = p^Co' = p,{dq' - ^'dt) (A.15) 

The restriction of the push-forward : T{A) T{ji{Vn+i)) to the vertical 
subspace V{A) determines a vector bundle homomorphism 



V{A) 
C 
A 



Composing with diagram ()A.7l a) and setting g := g ■ i^t, we get a further 
homomorphism 



V{h[Vn+l)) 

c 

Jl(Vn+l) 



V{A) 
C 

A 



(A.16) 



In coordinates, eqs. (IA.7b ). (IA.16|) provide the representation 



dz"^ 



\dz^)Xdq^M^^^ 



(A.17) 



Finally, by diagrams (fAl3]) . (fXSj) . setting k := K-i we get a bundle morphism 

^ l^*(V„+i) 



CM) 

c 

A 



(A.18) 



+1 



described in coordinates as 
t{k{a)) = t{a), 



i(k{a)) = q'{a), pi{k{a)) = pi{a) 



The latter allows to regard the contact bundle C{A) as a fibre bundle over the 
phase space l^*(V„+i), identical to the pull-back of y*(V„+i) through the map 

A^Vn+l. 
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B Finite deformations with fixed end points: an exis- 
tence theorem 



(i) Given an admissible, piecewise differentiable section 7 : [to, ii] ^ Vn+i, a crucial 
question is establishing under what circumstances every admissible infinitesimal 
deformation vanishing at the end points of 7 is tangent to an admissible finite 
deformation 7^ with fixed end points. The following preliminaries help simplifying 
the discussion. 

Lemma B.l Let 7 : (c, d) ^ A be the lift of an admissible differentiable section 
7 : (c, d) — > Vn+i ■ Then, for any closed interval [a, b] C (c, d) there exists a fibered 
local chart {U, k), k = {t, , . . . , q"^, z'^, . . . , z^) satisfying the properties 

j{t)eU yte[a,b]; 

the intersection 7((c, d)) D U coincides with the curve q^ = z^ = 0; 







V7(t) e U. 



(B.la) 
(B.lb) 
(B.lc) 



Proof. A straightforward argument, regarded as known, ensures the existence of 
local charts {V, h) , h = {t, f , . . . , g") in V„+i and {W, k') , k' = {t,x^ , . . . , 
in A satisfying the conditions 



7([a,6])cy, r(7(i)) =0 
7([a,6])cT^, x"(7(t)) = 



yt£ {c,d)n-f-\V) 

vt G {c,d)n^-\w) 



Without loss of generality we may assume 7r(VF) C V. The restriction to W 
of the projection vr : ^ — > Vn+i is then described in coordinates as 



with rank 



= n. In particular, the differentials dt , dq^, . . . , dq'^ are 

linearly independent everywhere on W. 

Let := ijL^a{t) dx'^\x^(t) denote r linear differential forms along 7, depending 
differentiably on t, and completing dg*|^(t) to a basis of T*^^^^{A) . 



Define r differentiable functions on W by 

n+r 



X^ 



a=l 



Then, by construction, the Jacobian 



is non singular at each point 



7(t) . The functions t, g*, z-^ form therefore a coordinate system in a neighborhood 
U of the intersection 7((c, d)) H W . The system is automatically fibered over 
Vn+i, and satisfies both properties (IB.ll a. b). and the first condition (IB. lb ) . 
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To complete the proof, let = 'ip^{t,q\z^) denote the representation of the 
imbedding A ji{Vn+i) in the coordinates t,q^,z^. Under an arbitrary linear 
transformation = a*j(t) , = z^ we have then the transformation laws 



In particular, if the matrix a^j{t) is a solution of the differential equation 

7(i) 



the coordinates t^q^,z^ satisfies all stated requirements. □ 

Every local chart ([/, A;) satisfying eqs. (jB.ll a. b, c) will be said to be adapted to 
the closed arc (7, [a, 6]) . 

Corollary B.l Let 7 = { ( 7'-''^ [a^-i, a^] ) , s = 1,...,A'"} be the lift of an ad- 
missible piecewise differentiable section (7, [to,ii])- Then, there exist fibered lo- 
cal charts {Us,ks), kg = [t, q^^y. . . ,q^y z^^y . . . , z^^-^) adapted to the arcs 7^*) 
such that, in each intersection tt{Us) n7r(?7c,+i), the coordinate transformation 
1{s+i) = 1{s) ' • • • ' ) satisfies the condition 



5) (B.2) 



The result follows at once from the proof of Lemma IB.ll . observing that the 
adapted coordinates are defined up to arbitrary linear transformations of the form 
qii _ j^i qj ^ z'^= z^, {A^- = const.), thus leaving full control on the values of 

the Jacobians ( — ^^i-^ at the corners. □ 



Every family |(C/s,/cs),s = 1,...,A^} satisfying the requirements of 
Corollarv IB . 1 1 will be said to be adapted to the lift 7. 

Assigning an adapted family of local charts automatically singles out a dis- 
tinguished infinitesimal control along each arc 7^**) , uniquely defined by the 
requirement 

In view of eqs. (j2.34l b). (|2.35l a). (|B.lb ). the absolute time derivative associated 
with /i(^) is described in coordinates as 



D d 



,N (B.3) 
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Noting that, by eq. (IB. 21) . the fields ( I are continuous at the corners 
7(0^) we conclude that the sections e(j) : [toj^i] — ^ ^(7) given by 

ewW=(7rrj yte[as-i,as], s = l,...,N (B.4) 

form a basis for the space Vh of /i-transported vector fields along 7. 

On account of eq. (jB.lb ). the corresponding dual basis for the space is given 
by eW(t) = = Vi G [0^-1,0^], s = l,...,Ar. By definition, 

together with eqs. ()B.3p we have therefore the dual relations 

(ii) Let us now come to the main question. Let 7 := { ( 7^^^ [a^-i, a^] ) } de- 
note an admissible, piecewise differentiable section, {(Us, kg)} a family of local 
charts adapted to 7, and {e(j)}, {e^*^} the corresponding dual bases for the spaces 

We recall that, with the notation of § 2.5, the most general infinitesimal de- 
formation X of 7 vanishing at t = to is determined by an element (Y, a) G 213 , 
namely by a vertical vector field Y along 7 and by a collection of real numbers 
a = (ai, . . . , aAT-i) . In particular, a necessary and sufficient condition for X 
to satisfy X{ti) = is expressed by the requirement ()2.53p which, in adapted 
coordinates, reads 

iy^{^)/^-i:-s[n7)U=0 (B.6) 

To inquire whether a given infinitesimal deformation vanishing at the end 
points of 7 is tangent to a finite deformation with fixed end points we intro- 
duce an auxiliary tool, namely a positive metric on Vh , described by a symmetric 
tensor $ = gij e^*) e^^^ . 

In view of the identification ^^(7) — [^Oj^i] ^ ^h, assigning $ automatically 
induces a scalar product along the fibres of V{^). This, in turn, determines a 
scalar product between vertical vector fields along 7, based on the identification 

{Y,Z) := {g{Y),m)) (B.7) 

Q ■ V{^) ^ V{'y) denoting the homomorphism (I2.16|) . In adapted coordinates, 
eqs. (iirril . (IRTI) provide the evaluation (Y, Z) = GabY^Z^, with 

As usual, the inverse of the matrix Gab will be denoted by G 




7/ ^ ^ 7 ^ ^7 
ab 
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By means of <i>, to every a = {ai, . . . ,ai\[-i) G ^ we associate — 1 
functions as(C) according to the prescription 

asiO ■■= as + asC - ^a,^^^9iji^'[i^'{l)],^ s = l,...,N-l (B.9) 
For notational convenience, the family is completed by the constant functions 

In a similar way, given any vertical vector field Y along 7, meant as a family 
of fields Y(s) = -(s) along the arcs of 7, for each € V/j we denote by 

fj^^'^^^ : tt{Us) Us, s = 1,...,N the (n + l)-parameter families of sections 
described in coordinates as 



It goes without saying that, being strictly coordinate-dependent, eq. ()B.10|) has 
no intrinsic geometrical meaning, but is merely a technical tool, whose usefulness 
will be clear in the subsequent discussion. 

On account of eqs. (jB.lOh . (IB.lip it is easily seen that, given any open subset 
A dVh with compact closure, there exists m > such that the image cTj.^*^^ (7r(f7s)) 
is entirely contained in Us for all 1^ ^ A, \^\ < m, s = 1, . . . ,N. 

Theorem B.l Let 7 be an admissible, piecewise differentiable evolution, and 
{Y, a) an admissible infinitesimal deformation of 7 vanishing at the endpoints. 
Define the metric <1> and the functions X(^yi{t) , o-siO ^■^ above. Then, given any 
open subset A C Vh with compact closure there exist an e > and a family 
^ii^^) ~ { ( ['^s-i(C); } of piecewise differentiable admissible sections 

defined for \^\ < e, G A and satisfying the properties 



(B.IO) 



with 




(B.ll) 



a) l{o,u)it) = lit) Vz^; 




the control equation 




(B.12) 
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Proof. Let A C Vh denote an open set with compact closure containing A. 

We choose m £ M4. as above, and examine the situation separately in each adapted 

chart {Us,ks). 

There, solving eq. (jB.12p amounts to determining the integral curves of the 
(n + l)-parameter family of vector fields ^^^^^^ ~ &i ~^ ^{s)'d(f '^(^s), with 

This, in turn, is equivalent to determining the integral curves of a single vector 
d d 

field Z^^-j = ^ + ^{s)'dq^ product manifold {—m,m) x j4 x tt{Us)- 

Let cl^^^^.^{t,x) denote the integral curve of Z(^) through the point {^,v,x). 
Also, let Xs-i denote the corner 7(as_i). Then, on account of eq. (|B.lb ). for any 
V* £ A the curve x^-i) coincides with the coordinate line = 0, ^ = 0, 

V = V* , and is therefore defined for all t in an open interval {bg-i-, hs) D [os-i, a^] . 

By well known theorems in ordinary differential equations [13^ [3] this im- 
plies the existence of an open neighborhood Wg-i 9 (0,i^*, x^-i) such that the 
curve cl^^^^^{t,x) is defined for all (^,1/, x) G Ws-i and all t in the closed interval 

[t{x),as{0] C {bs-i,bs). 

In particular, denoting by the slice t = as(^) in (—m,m) x ^ x tt{Us), we 
conclude that the 1-parameter group of diffeomorphisms determined by the field 
Zf^g-j maps the intersection Wg-i H T,s-i into an open neighborhood of the point 
(0, u*, Xs) in . Without loss of generality we may always arrange that the image 
of each Ws-^i Ci S^-i is contained in Ws H , s = 1, . . . , N . 

The rest is now entirely straightforward: let U and e^; > respectively 
denote an open neighborhood of u* in A and a positive number such that 
(^,z/, xo) € WoHHo V 1^1 < , € [/ (notice that, according to our thesis, we are 
"freezing" the choice of the point xq). For each |^| < Eu , f € U consider the 
sequence of closed arcs '■ [q^s-i(0> ~^ t^{Us) defined inductively by 

The collection 7(^_t,) := {{"f^^^^y[as-i{0,as{0])^ « = l,---,^} is then easily 
recognized to define an (n + l)-parameter family of continuous, piecewise differ- 
entiable sections satisfying all requirements a), b), c), d) of the Theorem. 

To complete our proof let us finally recall that, for any u* £ A, the family 7(g,jy) 
exists for all i' in an open neighborhood U 3 v* and all |^| < Eu . On the other 
hand, by the assumed compactness of A, the subset A C A may be covered by a 
finite number of subsets {[/i, . . . , t/fc} of the required type. 

This establishes Theorem IB. I} with e = min {e^^ , • • • , Eu^ } . □ 
According to Theorem IB. 11 for any open subset A C V/i with compact clo- 
sure, the correspondence v — > 7(g,;y)(ti) sets up a 1-parameter family of differen- 
tiable maps of A into the hypersurface t = ti, with values in a neighborhood 
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of the point 7(^1)- Moreover, given any differentiable curve v = v{C) in A, 
the 1-parameter family of sections 7(^,1/(5)) (*) > I'^l < e, i G [^Oi^i] is a defor- 
mation of 7 tangent to the original infinitesimal deformation X determined by 
(y, aN-i), and leaving the first end point xq fixed. 

Our original problem is therefore reduced to establishing the existence of a 
curve satisfying 7(g = 7(ti) in some open neighborhood of ^ = 0. 

In adapted coordinates, setting for simplicity := V'(7v)(^i' ^' ^) ' 

required condition reads 

= i = l,...,n (B.13a) 

Taking the relations V9*(0, z^) = g()^)(7(ti)) = 0, (^)^^q = X^ti) into ac- 
count, a straightforward application of Taylor's theorem shows that, whenever the 
condition X{ti) = holds true, i.e. whenever the field Y and the coefficients Ug 
satisfy eq. (|B.6p . the functions y?* are necessarily of the form (/?*(^, i/) = i^) , 

with 6'*(^, v) regular at ^ = 0. Under the stated assumptions, eq. (|B.13l a) is there- 
fore equivalent to the condition 

0*(^,zy\...,i/") = i = l,...,n (B.13b) 

Let us discuss the solvability of eqs. (|B.13l b) for the z/*'s as functions of 
in a neighborhood of ^ = 0. To this end with we observe that the matching 
conditions c) of Theorem 13.11 give rise to relations of the form 

^{s+i) ~ ^(s+i)(^' 1{s) ' • • • ! denoting the transformation between adapted co- 
ordinates in the intersection ti{Us H [7^+1 )• From these, deriving with respect to ^ 
we get the expressions 



dt di di dt di dq^, V dt di 

At ^ = 0, recalhng eqs. (plil]) . (|R2]) . ([R9]) as well as the identification X 



the latter provide the relation 



dq 



(s+i) 



(B.15) 



In a similar way, on account of eqs. (IB.2p . (|B.9p . (IB.15|) . deriving eq. ()B.14h 
with respect to ^ and evaluating everything at ^ = , a straightforward calculation 
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yields the result 



2^ i 



a. 



(s+i) X ^ X 



2 a, 



expressing the jumps 



dt dt 



+ a.. 



grk 



V''^(7) u'' (B.16) 



in terms of the section 7, of the 



infinitesimal deformation and of the variables . 

In addition to this let us now make use of the fact that, in each adapted chart, 
eqs. (|B.12p imply the evolution equations 



+ 



dq^dq^ 



X'^X'' + 2 



■As) 



Q2^i 



QqkQ^A 



y(s) 




the cancelation arising from eq. (|B.lb ). 

Prom the latter, restoring the notation ip^{^,v) for ipi^j^-^{ti,^,v) and recalling 
eqs. (jB.lip . ()B.16p . as well as the fact that the components gij are by definition 
constant along 7 , we get an expression of the form 



as / Q^i 



N-l 



9jk 



^'{1) Igjki^" (B.17) 



with 6' G M depending only on the section 7 and on the original infinitesimal 
deformation. Collecting all results we can therefore state 

Proposition B.l Let 7 : [toi^i] ~^ V^+i be a continuous, piecewise differentiahle, 
admissible section. Then, if the matrix 



to 



\dz^)Adz'^ 



N-l 



dt + Y,^s V'^'(7) 



(B.18) 



s=l 
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is non-singular, every infinitesimal deformation of 7 vanishing at the end points 
is tangent to a finite deformation with fixed end points. 

Proof. The conclusion follows at once from the fact that, on account of eq. (IB.17|) . 
the non-singularity of the matrix (jB.lSp ensures the solvability of eqs. ()B.13b ) in 
a neighborhood of ^ = . □ 

Proposition lB.ll mav be rephrased in the language of § 2.5 observing that, whenever 
the section 7 is abnormal, the matrix (|B.18|) is necessarily singular. 

Under the stated circumstance, in fact. Proposition 12.51 and eq. ()B.4p imply the 
existence of at least one non-zero virtual 1-form piCb^\^ with constant components 
Pi obeying the relations 

and therefore automatically satisfying pi S^^ = . 

More specifically, denoting by p the abnormality index of 7, we have the fol- 
lowing 

Theorem B.2 The matrix iB.18\) has rank n — p. 

Proof. By definition, the index p coincides with the dimension of the annihilator 
(T(2IJ))'^ C , which, in turn, is identical to the dimension of the space of 
constant solutions of eqs. (|B.19p . 

On the other hand, by eqs. (jB.Sh . (IB.lSp . the matrix S^^ is positive semidefinite. 
Its kernel is therefore identical to the totality of zeroes of the quadratic form 
S^-' PiPj , i.e. to the totality of n-tuples (pi , . . . , Pn) ^ satisfying the relation 

Due to the positive definiteness of G'^^{t), the last condition is equivalent to 
eqs. (|B.19p . This proves dim ( ker(5*-')) = p ^ rank(S'*-') =n — p. □ 

In the language of § 2.4, Proposition IB . 1 1 and Theorem lB . 21 show that the normal 
evolutions form a subset of the ordinary ones, thus establishing Proposition 12.61 
Along the same lines, a deeper result is provided by the following 

Theorem B.3 Let p (> 0) denote the abnormality index of ^ . Then a sufficient 
condition for every infinitesimal deformation vanishing at the end points of 7 to be 
tangent to a finite deformation with fixed end points is the existence of an (n—p)- 
dimensional submanifold S C V^+i contained in the slice t = ti and containing the 
point 7(ti), such that every deformation 7^ leaving 7(to) fixed satisfies 7g(ti) S S 
for all sufficiently small. 
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Proof. Assuming the existence of a submanifold S — > V„+i with the stated 
properties, we denote by {V,(^,---, C"^"") a local chart in S centered at the point 
7(ti), and by 

t = h, q^l,) = g'{C\...:C-'') (B.20) 

the representation of S in adapted coordinates. 

Then, for any open subset A G Vh with compact closure and for sufficiently 
small ^, the correspondence — > 7(5,i/)(ii) factors through S, giving rise to a 
differentiable map g : (— e, e) x A ^ S satisfying the relation 7(^^,y)(ti) = i- giS,-, f)- 

In coordinates, setting C°'{g{S,, i^)) = i^^, ■ ■ ■ i and resuming the nota- 
tion i/^, . . . , z^") for '?(-]y)(7(5.j^)(*i) 5 this provides the identification 

= Q\g\i,v\...,v^),...,g^~^{i,u\...,v^)) (B.21) 

Prom the latter, recalling the relation (^"(O, z^"*^, . . . , z/") = C"(7(^i)) = as 

well as the fact that the Jacobian g^^i ...^n-p) has maximal rank it is easily seen 

that the equalities ip^{0,i'^, ... ,!/"■) = ^(0, u^, . . . , v^) = are reflected into 
analogous properties of the functions g°' . 

By Taylor's theorem we have therefore an expression of the form 

5° = eV"(e,i^\...,i^") (B.22) 

with the functions /i" regular at ^ = 0. 

The proof is thus reduced to establishing the solvability of the system 

/x°(e,zv\...,z^") = (B.23) 

for the i^*'s as functions of ^ in a neighborhood of ^ = 0. 

To this end, by direct computation, from eqs. ()B.2ip . (|B.22p we derive the 
relation 

Together with eqs. ()B.17p . (|B.18p . the latter provides the identification 

b^ + S^'grk^^' = 2(^) /x-(0,z.\...,z.") (B.24) 

In view of eq. dEMl), the functions /i°(0, z/^, . . . , i^") are linear polynomials 

^"(0, = M"fc u'' + c" (B.25) 

with coefficients M°^k, uniquely determined in terms of 6*, S^^ , grk and of the 
imbedding (|B.20p . In particular, by eq. ()B.24p . the rank of the matrix M^^ cannot 
be smaller than the rank of S^^ and, of course, cannot exceed n — p. According to 
Theorem IB.2i we have therefore rankM^^ = n — p. 
Collecting all results we conclude: 
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the system ()B.23P admits oo^' solutions of the form (0, v*^ , . . . , i 



on account of eq. (|B.25p . the Jacobian 



each point (0, , . 



has rank n — p at 



By continuity, it has therefore rank n — p in a 



neighborhood of every solution (0, v 



*i 



of eqs. (Ib:23]) . 



By the implicit function theorem, this proves that the system ()B.23p admits at 
least a solution of the form = in a neighborhood of ^ = (actually, 

infinitely many solutions whenever p > Q). □ 
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